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I. INTRODUCTION 



We discuss the possible relevance of gravitational-wave (GW) experiments for physics at very 
high energy. We examine whether, from the experience gained with the computations of various 
specific relic GW backgrounds, we can extract statements and order of magnitude estimates that are 
as much as possible model-independent, and we try to distinguish between general conclusions and 
results related to specific cosmological mechanisms. We examine the statement that the Virgo/LIGO 
experiments probe the Universe at temperatures T ~ 10 7 — 10 10 GeV (or timescales t ~ 10 -20 — 10~ 26 
sec) and we consider the possibility that they could actually probe the Universe at much higher 
energy scales, including the typical scales of grand unification, string theory and quantum gravity. 
We consider possible scenarios, depending on how the inflationary paradigm is implemented. We 
discuss the prospects for detection with present and planned experiments. In particular, a second 
Virgo interferometer correlated with the planned one, and located within a few tens of kilometers 
from the first, could reach an interesting sensitivity for stochastic GWs of cosmological origin. 
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The energy range between the grand unification scale Mqtjt ~ 10 16 GeV and the Planck scale Mpi ~ 1.22 x 10 19 
GeV is crucial for fundamental physical questions and for testing current ideas about grand unification, quantum 
gravity, string theory. Experimental results in this energy range are of course very difficult to obtain. From the 
particle physics point of view, there are basically two important experimental results that can be translated into 
statements about this energy range (see e.g. jl]] for recent reviews): (i) the accurate measurement of gauge coupling 
Vh constants at LEP that, combined with their running with energy, shows the unification of the couplings at the scale 
jpP" A/gut, provided that the running is computed including the supersymmetric particles in the low energy spectrum. 
J> . And (ii) the negative results on proton decay; the lower limit on the inverse of the partial decay width for the processes 
Is^J ' p — ► e + 7r° and p — > K + v are 5.5 x 10 32 yr and 1.0 x 10 32 yr, respectively, and imply a lower bound Mqut^IO 15 GeV, 
which excludes non-supersymmetric SU(5) unification. Further improvement is expected from the SuperKamiokande 
experiment, which should probe lifetimes ~ 10 34 yr. 

From the cosmological point of view, informations on this energy range can only come from particles which decoupled 
from the primordial plasma at very early time. Particles which stay in thermal equilibrium down to a decoupling 
temperature Td oc can only carry informations on the state of the Universe at E ~ Idee- ^ u informations on physics 
at higher energies has in fact been obliterated by the successive interactions. 

The condition for thermal equilibrium is that the rate T of the processes that mantain equilibrium be larger than 
the rate of expansion of the Universe, as measured by the Hubble parameter H ||. The rate is given by T = na\v\ 
where n is the number density of the particle in question, and for massless or light particles in equilibrium at a 
temperature T, n ~ T 3 ; \v\ ~ 1 is the typical velocity and a is the cross-section of the process. Consider for instance 
the weakly interacting neutrinos. In this case the equilibrium is mantained, e.g., by electron-neutrino scattering, and 
at energies below the W mass a ~ G 2 F (E 2 ) ~ G 2 F T 2 where Gf is the Fermi constant and (E 2 ) is the average energy 
squared. The Hubble parameter during the radiation dominated era is related to the temperature by H ~ T 2 /Mp\. 
Therefore 



T ^ G 2 F T 5 ^ ( T \ 



H J ncutlino ~ T 2 /M P r \1McV J 



(1) 



Even the weakly interacting neutrinos, therefore, cannot carry informations on the state of the Universe at temper- 
atures larger than approximately 1 MeV. If we repeat the above computation for gravitons, the Fermi constant Gf 
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is replaced by Newton constant G = 1/Mp, (we always use units H = c = ks = 1) and at energies below the Planck 
mass 



The gravitons are therefore decoupled below the Planck scale. (At the Planck scale the above estimate of the cross 
section is not valid and nothing can be said without a quantum theory of gravity). It follows that relic gravitational 
waves are a potential source of informations on very high-energy physics. Gravitational waves produced in the very 
early Universe have not lost memory of the conditions in which they have been produced, as it happened to all 
other particles, but still retain in their spectrum, typical frequency and intensity, important informations on the 
state of the very early Universe, and therefore on physics at correspondingly high energies, which cannot be accessed 
experimentally in any other way. It is also clear that the property of gravitational waves that makes them so interesting, 
i.e. their extremely small cross section, is also responsible for the difficulties of the experimental detection.^] 

With the very limited experimental informations that we have on the very high energy region, Mqjjt<E<Mp\, it 
is unlikely that theorists will be able to foresee all the interesting sources of relic stochastic background, let alone to 
compute their spectra. This is particularly clear in the Planckian or string theory domain where, even if we succeed 
in predicting some interesting physical effects, in general we cannot compute them reliably. So, despite the large 
efforts that have been devoted to understanding possible sources, it is still quite possible that, if a relic background of 
gravitational waves will be detected, its physical origin will be a surprise. In this case a model-independent analysis 
of what we can expect might be useful. 

In this paper we discuss whether, from the experience gained with various specific computations of relic backgrounds, 
it is possible to extract statements or order of magnitude estimates which are as much as possible model- independent. 
These estimates would constitute a sort of minimal set of naive expectations, that could give some orientation, 
independently of the uncertainties and intricacies of the specific cosmological models. We discuss typical values of the 
frequencies involved and of the expected intensity of the background gravitational radiation, and we try to distinguish 
between statements that are relatively model-independent and results specific to given models. 

The paper is written having in mind a reader interested in gravitational-wave detection but not necessarily com- 
petent in early Universe cosmology nor in physics at the string or Planck scale, and a number of more technical 
remarks are relegated in footnotes and in an appendix. We have also tried to be self-contained and we have attempted 
to summarize and occasionally clean up many formulas and numerical estimates appearing in the literature. The 
organization of the paper is as follows. In sect. 2 we introduce the variables most commonly used to describe a 
stochastic background of gravitational waves. We give a detailed derivation of the relation between exact formulas 
for the signal-to-noise ratio, and approximate but simpler characterizations of the characteristic amplitude and of 
the noise. The former variables are convenient in theoretical computations while the latter arc commonly used by 
experimentalists, so it is worthwhile to understand in some details their relations. In sect. 3 we apply these formulas 
to compute the sensitivity to a stochastic background that could be obtained with a second Virgo interferometer 
correlated with the first, and we compare with various others detectors. We find that in the Virgo- Virgo case the 
noise which would give the dominant limitation to the measurement of a stochastic background is the mirror thermal 
noise, and we give the sensitivity for different forms of the relic GW spectrum. In sect. 4 and 5 we discuss estimates 
of the typical frequency scales. We examine the statements leading to the conclusion that Virgo/LIGO will explore 
the Universe at temperatures T ~ 10' GeV (sect. 4). and in the appendix we discuss some qualifications to this 
statement. In sect. 5 we discuss the possibility to reach much higher energy scales, including the typical scales of 
grand unification and quantum gravity. In sect. 6 we discuss different scenarios, depending on how the inflationary 
paradigm is implemented. Characteristic values of the intensity of the spectrum and existing limits and predictions 
are discussed in sect. 7. Sect. 8 contains the conclusions. 



'Thinking in terms of cross-sections, one is lead to ask how comes that gravitons could be detectable altogheter, since the 
graviton-matter cross section is smaller than the neutrino-matter cross section, at energies below the VK-mass, by a factor 
G 2 /G% ~ 10 -67 and neutrinos are already so difficult to detect. The answer is that gravitons are bosons, and therefore their 
occupation number per cell of phase space can be nt > 1; we will see below that in interesting cases, in the relic stochastic 
background we can have nk ~ 10 40 or larger, and the squared amplitude for exciting a given mode ot the detector grows as n\. 
So, we will never really detect gravitons, but rather classical gravitational waves. Neutrinos, in contrast, are fermions and for 
them rife < 1. 




(2) 
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II. DEFINITIONS 



A. n gw (/) and the optimal SNR 

The intensity of a stochastic background of gravitational waves (GWs) can be characterized by the dimensionless 
quantity 

V s Af) = -^- f , (3) 
Pc a log/ 

where p gw is the energy density of the stochastic background of gravitational waves, / is the frequency (lo = 2nf) 
and p c is the present value of the critical energy density for closing the Universe. In terms of the present value of the 
Hubble constant H$, the critical density is given by 

(4) 

The value of Hq is usually written as Hq = ho x 100 km/(sec-Mpc), where ho parametrizes the existing experimental 
uncertainty. Ref. || gives a value 0.5 < h < 0.85. In the last few years there has been a constant trend toward lower 
values of h and typical estimates are now in the range 0.55 < h < 0.60 or, more conservatively, 0.50 < h < 0.65. 
For instance ref. Q], using the method of type IA supernovae, gives two slightly different estimates h = 0.56 ± 0.04 
and h = 0.58 ± 0.04. Ref. f|, with the same method, finds ho — 0.60 ± 0.05 and ref. H, using a gravitational lens, 
finds ho — 0.51 ± 0.14. The spread of values obtained gives an idea of the systematic errors involved. 

It is not very convenient to normalize /? gw to a quantity, p c , which is uncertain: this uncertainty would appear in 
all the subsequent formulas, although it has nothing to do with the uncertainties on the GW background. Therefore, 
we rather characterize the stochastic GW background with the quantity /i§f2 gw (/), which is independent of ho- All 
theoretical computations of a relic GW spectrum are actually computations of dp gw /dlogf and are independent of 
the uncertainty on Hq. Therefore the result of these computations is expressed in terms of ft,Qf2 gw , rather than of 
il gw -0 

To detect a stochastic GW background the optimal strategy consists in performing a correlation between two (or 
more) detectors, since, as we will discuss below, the signal will be far too low to exceed the noise level in any existing 
or planned single detector (with the exception of the space interferometer LISA, see below). The strategy has been 
discussed in refs. [|7[— |l0|] , and a clear review is ref. jy]. Let us recall the main points of the analysis. The output of 
any single detector is of the form s;(t) — hi(t) + rii(t), where i = 1, 2 labels the detector, and the output Si(t) is made 
up of a noise rii(t) and possibly a signal hi(t). In the typical situation, hi <C n,;. We can correlate the two outputs 
defining 

S= [ dt [ dt' s 1 {t)s 2 {t')Q{t - 1') , (5) 

J-T/2 J-T/2 

where T is the total integration time (e.g. one year) and Q a filter function. If the noises in the two detectors are 
uncorrelated, the ensemble average of the Fourier components of the noise satisfies 

(n*(f)n j (f))=S(f-f')S ij ^H\f\). (6) 

The above equation defines the functions Sn\\f\), with dimensions Hz -1 . The factor 1/2 is conventionally inserted 
in the definition so that the total noise power is obtained integrating S n (f) over the physical range < / < oo, 

rather than from — oo to oo. The noise level of the detector labelled by i is therefore measured by h$ = y Sn ^ , with 
dimensions Hz -1 / 2 . The function S n is known as the square spectral noise density.^] 



^This simple point has occasionally been missed in the literature, where one can find the statement that, for small values of 
Ho, f2 gw is larger and therefore easier to detect. Of course, it is larger only because it has been normalized using a smaller 
quantity. 

* Unfortunately there is not much agreement about notations in the literature. The spectral noise density, that we denote by 
S n (f) following e.g. ref. §, is called P(f) in ref. [0. Other authors use the notation Sh(f), which we instead reserve for the 
spectral density of the signal. To make things worse, S„ is sometime defined with or without the factor 1/2 in eq. (0). 
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As discussed in refs. |TT[, in the limit hi <C rii, for any given form of the signal, i.e. for any given functional form 
of /&o^gw(/)) it is possible to find explicitly the filter function Q(t) which maximizes the signal-to-noise ratio (SNR). 
In the case of L-shaped interferometers the corresponding value of the optimal SNR turns out to be (see e.g. ref. [ pj , 
eq.(43)) 



SNR 



T 



df 



7 2 (/)^ w 



(/) 



1/4 



(7) 



We have taken into account the fact that what has been called S in eq. (^|) is quadratic in the signals and, with usual 
definitions, it contributes to the SNR squared. This differs from the convention used in ref. The function j(f) 
is called the overlap function. It takes into account the difference in location and orientation of the two detectors. It 
has been computed for the various pairs of LIGOl, LIG02, Virgo and GEO detectors ||. For detectors very close 
and parallel, 7(/) = 1. Basically, j(f) cuts off the integrand in eq. (0) at a frequency 2nf of the order of the inverse 
separation between the two detectors. For the two LIGO detectors, this cutoff is around 60 Hz. We will discuss 7(/) 
in sect. 2C, where we will also comment on the modifications needed for different geometries. 

In principle the expression for the SNR, eq. (]?]), is all that we need in order to discuss the possibility of detection 
of a given GW background. However it is useful, for order of magnitude estimates and for intuitive understanding, 
to express the SNR in terms of a characteristic amplitude of the stochastic GW background and of a characteristic 
noise level, although, as we will see, the latter is a quantity that describes the noise only approximately, in contrast 
to eq. (^) which is exact. We will introduce these quantities in the next two subsections. 



B. The characteristic amplitude 

A stochastic GW at a given point x = can be expanded, in the transverse traceless gauge, as (we follow the 
notations of ref. plj ], app.A) 

h ab (t) = ]T f df f dClh A (f,Cl)exp(2mft)e±(n), (8) 

where /, fi) = h A (f,Cl). 17 is a unit vector representing the direction of propagation of the wave and dfl = 
dcos0d(f>. The polarization tensors can be written as el 6 (fi) = m a mb — n a ni, and e* 6 (f2) = m tt n& + n a m , with m, n 
unit vectors ortogonal to f2 and to each other. With these definitions, e A j JQ)e A ' ab (tl) = 25 AA . For a stochastic 
background, assumed to be isotropic, unpolarized and stationary (see |jn]jL2| for a discussion of these assumptions) 
the ensemble average of the Fourier amplitudes can be written as 

(^(/,o)^(^nO) = ^(/-/0^ 2 (^^ , )W^ fc (/) ) (9) 

where 5 2 (tt, O') = 8{<f) — (j)')8{cos6 — cos 9'). The function £/,(/) defined by the above equation has dimensions Hz -1 
and satisfies Sh(f) = Sh(—f)- The factor 1/2 is conventionally inserted in the definition of Sh in order to compensate 
for the fact that the integration variable / in eq. (§h ranges between — oo and +oo rather than over the physical 
domain < / < oo. The factor l/(47r) is inserted so that, integrating the left-hand side over dCt and over dCt' , we 
get S(f — f')5AA'(l/2)Sh(f)- With this normalization, Sh(f) is therefore the quantity to be compared with the noise 
level S n (f) defined in eq. (g). Using eqs. (||J|) we get 

/OO /•/ — oo 

dfS h (f) = 4 d(log/) fS h (f). (10) 

-oo Jf=0 

We now define the characteristic amplitude h c (f) from 

rf=oo 

(h ab (t)h ab (t))=2 d(\ogf)h 2 c (f). (11) 

Jf=0 

Note that h c (f) is dimensionless, and represents a characteristic value of the amplitude, per unit logarithmic interval of 
frequency. The factor of two on the right-hand side of eq. (pd[) is inserted for the following reason. The response of the 
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detector to a single wave with amplitudes h + ,h x is of the form (Ref. jl3|, eqs. (26) and (103)) h(t) — F + h + +F x h x . 
For an interferometer, h(t) = Al(t)/L where I is the difference in arm lengths. The function h(t) is also called the 
gravitational wave strain acting on the detector. The functions f+, x are known as detector pattern functions, and 
< l-F+.x I < 1- They depend on three angles, that determine the direction of arrival of the GW and its polarization. 
For any quadrupole-beam-pattern GW detector we have (ref. O], pg. 369) (F 2 )^ = (F x )q an d (F+Fx)q = 0, where 

V ln J 4tt 2tt V ' y ' 

denotes an average over the direction of propagation of the wave, dCt = d cos 6d(f>, and over the angle ip that gives 
the preferred frame (x',y') where the wave in the transverse traceless gauge takes the simple form h x 'x' = ~h y l y l , 
h x t y t = h y t x t (see jl3]], pg. 367). This average should not to be confused with (. . .) which is the time average, i.e., in 
Fourier space, the ensemble average of eq. (^). Therefore, {h 2 {t))^ = {F] r )^h 2 + + {F x )^ l h 2 x . In a unpolarized stochastic 
background (h 2 ^) = {h 2 x ) (where (. . .) is the ensemble average of eq. (^)) and therefore, since (F 2 )^ — (F 2 )^, the two 
terms in h 2 (t), after averaging over the angles and over time, give the same contribution. This motivates the factor 
of two in the definition of h 2 (f), eq. (|ll|). For the same reason, we could insert a factor (F 2 )^ in the definition of 
h 2 (f). However, what is really meaningful is not a characterization of the signal h c (f) nor of the noise level h n (f) 
separately, but rather the signal-to-noise ratio discussed in the previous subsection. We are trying to express the SNR 
in terms of a ratio of two quantities, h c (f)/h n (f)- So, we are free to arbitrarily move factors from h c to h n as long as 
h c /h n is unchanged. However, it is convenient to collect in h c all factors related to the source, and in h n all factors 
related to the detectors. We do not insert a factor (F 2 ) 1 ^ 2 in the definition of h c , and therefore, automatically, we 
will obtain a factor (-F+)^ in h n , see eqs. ( [23|]2"e| ) below. The same is true for numerical factors like the factor of 

two just discussed. We could as well have neglected it in h 2 and we would find an additional factor l/v2 in h n . We 
will discuss h n in the next subsection. 
Comparing eqs. ( |l0| ) and Q), we get 

h 2 c (f) = 2fS h {f). (13) 

We now wish to relate h c (f) and /iQfi gw (/). The starting point is the expression for the energy density of gravitational 
waves, given by the 00-component of the energy-momentum tensor. The energy-momentum tensor of a GW cannot 
be localized inside a single wavelength (see e.g. ref. Il5|, sects. 20.4 and 35.7 for a careful discussion) but it can be 
defined with a spatial averaging over several wavelenghts: 

' (h ab h ab ) . (14) 



32ttG 



For a stochastic background, the spatial average over a few wavelengths is the same as a time average at a given point, 
which, in Fourier space, is the ensemble average performed using eq. (^]). We therefore insert eq. (^) into eq. ( |l4| ) and 
use eq. (0). The result is 



32ttG 

so that 



/— oc 

Pgw = tjttt; / d ( lo §/) f(^f) 2 S h (f) , (15) 
f=o 



dlogf 2G 



7^fS h {f) • (16) 



Comparing eqs. fllq) and (13) we get the important relation 



= ^fhl{f), (17) 



or, dividing by the critical density p c , 



dlogf AG 



^gw(/) = |^/ 2 ^(/). (18) 



Inserting the numerical value of H , we find (ref. EJ, eq. (65)) 
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h c (f) ~ 1.263 x 10~ 18 



(19) 



Using eqs. ([L3|,|l8]) we can also write f2 gw (/) = (An 2 /3H 2 )f 3 Sh{f)- Using this relation, and defining S n (f) 
(Sn\f)Sn\f)) 1 / 2 , eq. (0) can be written in a more transparent form, 



SNR 



. >TI df (wiy?m 



1/4 



(20) 



The physical reason for the appearance of the factor 2/5 in the above formula will be clear from eq. ( 53 ) below. (The 
number 2/5 is specific to L-shaped interferometers, see below). The factor of 2 in front of the integral can instead be 
understood from J_ oa df = 2 J Q df. 

Finally, we mention another useful formula which expresses /io^gw(/) m terms of the number of gravitons per 
cell of the phase space, n(x,k). For an isotropic stochastic background n(x, k) = rif depends only on the frequency 
/ = |fc|/(27r), and p gw = 2 / n f 2wf d 3 k/(2n) 3 = 16tt 2 / °° d(log f)n s f. Therefore dp gw /d\ogf = IGtt 2 ^/ 4 , and 



,(/)^3.e(- 



10 37 / V 1kHz 



(21) 



As we will discuss below, to be observable at the LIGO/Virgo interferometers, we should have at least /i 2 ,^ 
between 1 Hz and 1 kHz, corresponding to n/ of order 10 31 at 1 kHz and nt ~ 111 1 ' 
GW background is therefore exceedingly classical, rik ^> 1. 
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10 43 at 1 Hz. A detectable stochastic 



C. The characteristic noise level 



We have seen in the previous section that there is a very natural definition of the characteristic amplitude of the 
signal, given by h c (f), which contains all the informations on the physical effects, and is independent of the apparatus. 
We can therefore associate to h c (f) a corresponding noise amplitude h n (f), that embodies all the informations on the 
apparatus, defining h c (f)/h n (f) in terms of the optimal SNR. 

If, in the integral giving the optimal SNR, eq. (0) or eq. (p0[), we consider only a range of frequencies A/ such that 
the integrand is approximately constant, we can write 



SNR: 



8 ^a,7 2 (/)^(/) 



—TAJ 
25 



S 2 n (f) 



1/4 



2TA/ 7 2 (/)^(/) 
25/ 2 S2(/) 



1/4 



(22) 



The right-hand side of eq. ( |22| ) is proportional to h c (f), and we can therefore define h n (f) equating the right-hand 
side of eq. (p2|) to h c (f)/h n (f), so that 



M/) 



(±TA/)V4 



-JSn(f) 



1/2 



1(f) 



1/2 



(23) 



For L-shaped interferometers, the overlap function j(f) is defined in terms of the detector pattern functions F+\, 
where i — 1, 2 labels the detector, as M 



-in 



7Cf)=5/^«* rt/n - A *W ) 4 a) +^: 



.(1) F (2) 



(24) 



If the separation Ax between the two detectors is very small, 2irfAx <C 1 and if the detectors have the same 
orientation, so that F^ x — -Fl x then (averaging also over the angle ip discussed in the previous subsection) 



(25) 



For L-shaped interferometers (F 2y 



- 1/5 (see ref. |J, eq. (110), or ref. Q). We see that the factor 5/2 in the 
definition of j(f), eq. (p3), has been inserted so that for parallel detectors at the same site ^(f ) = 1. Therefore, the 
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factor 5 1 / 2 in eq. (|2^) measures the increase in the noise level due to the fact that stochastic GWs hit the detectors 
from all directions, rather than just from the direction where the sensitivity is optimal, while I/7 1 / 2 measures the 
decrease in sensitivity due to the detectors separation. 

For a generic geometry, the factors 5 in the above formulas are replaced by l/{Fl)^. In the special case j(f ) = 1 

(detectors at the same site) the factor (5/7) 1 / 2 in eq. ((2^) becomes therefore 1/(F 2 _)V 2 , and we recover eq. (66) of 
ref. |l|] (for comparison, note that the quantity denoted Sh(f) in jl3j is called here S n (f)/2). 

From the derivation of eq. ( p3[ ) we can understand the limitations implicit in the use of h n (f). It gives a measure 
of the noise level only under the approximation that leads from eq. (pp|), which is exact (in the limit hi <C rii), to 
eq. (|22|). This means that A/ must be small enough compared to the scale on which the integrand in eq. ( J20| ) changes, 
so that r y(f)Sh(f)/S n (f) is approximately constant. In a large bandwidth this is non trivial, and of course depends 
also on the form of the signal; for instance, if ft.Q^ g w is flat, then Sh(f) ~ l// 3 - For accurate estimates of the SNR 



at a wideband detector there is no substitute for a numerical integration of eq. (f7|) or eq. (20). However, for order 



of magnitude estimates, eq. (19) for h c (f) and eq. ( p3| ) for h n (f) are simpler to use, and they have the advantage 
of clearly separating the physical effect, which is described by h c (f), from the properties of the detectors, that enter 
only in h n (f). 

Eq. (|2^) also shows very clearly the advantage of correlating two detectors compared with the use of a single 
detector. With a single detector, the minimum observable signal, at SNR=1, is given by the condition Sh(f) > S n (f). 
This means, from eq. (fi~3|), a minimum detectable value for h c (f) given by h^ in = (2/S , n (/)) 1//2 . The superscript Id 
reminds that this quantity refers to a single detector. From eq. ( |23| ) we find for the minimum detectable value with 
two interferometers in coincidence, /i 2 f in , 

1/2 / 1 tt \ 1/4 /, s 1/4 -1 



C in (/) = ? T ? ^^(/)(^) - 1-77 x lO-XK/) (^) ^y- (26) 

Of course, the reduction factor in the noise level is larger if the integration time is larger, and if we increase the 
bandwidth A/ over which a useful coincidence (i.e. 7(/) ~ 1) is possible. Note that /i 2 f2 gw is quadratic in h c (f), so 
that an improvement in sensitivity by two orders of magnitudes in h c means four orders of magnitude in h^fl^. 



III. APPLICATION TO VARIOUS DETECTORS 



Single detectors. To better appreciate the importance of correlating two detectors, it is instructive to consider 
first the sensitivity that can be obtained using only one detector. In this case a hypothetical signal would manifest 
itself as an eccess noise, and should therefore satisfy Sh(f)^S n (f). Using eqs. ( |l3||llj| ) we write Sh in terms of /igf2 gw , 
and we write S n (f) = h 2 . For the minimum detectable value of h\fl m we get 

h « n ^ f ^ ir ims)\w^f ■ ,27 » 

Fig. 1, taken from ref. p6[ |, shows the sensitivity of the planned Virgo interferometer, i.e., the quantity hf as a function 
of /. We see that Virgo, used as a single detector, can reach a minimum detectable value for /igOgw of order 1CT 2 or 
at most a few times 10 -3 , at / = 100Hz. Unfortunately, this is not an interesting sensitivity level; as we will discuss 
in sect. 7, an interesting sensitivity level for /ip^gw is at least of order — 10~ 6 . To reach such a level with a single 
Virgo interferometer we need, e.g., h f < 3 x 1(T 25 - 10" 24 Hz" 1/2 at / = 100 Hz, orh f < 10~ 26 - 3 x 10" 26 Hz _1/2 
at / = 1 kHz. We see from fig. 1 that such small values of hf are very far from the sensitivity of first generation 
interferometers, and are in fact even well below the limitation due to quantum noise. 

If we consider instead the NAUTILUS detector, which has a target sensitivity hf = 8.6 x 10~ 23 Hz _1 / 2 and operates 
at two frequencies in the kHz region (907 Hz and 923.32 Hz) ]l7[] , we find ft-o^gw' ~ 5. The analysis using existing 
EXPLORER data (from 1991 and from 1994) and more recent NAUTILUS data gives a bound 0(300) on ft,^ gw @. 

Thus, we see that in these cases we cannot get a significant bound on /i 2 r2 gw . A very interesting sensitivity, possibly 
even of order /igfigw ~ 10~ 13 , could instead be reached with a single detector, with the planned space interferometer 
LISA @, at / - 10- 3 Hz. 

Note also that, while correlating two detectors the SNR improves with integration time, see eq. (Q), this is not so 
with a single detector. So, independently of the low sensitivity, with a single detector it is conceptually impossible 
to tell whether an eccess noise is due to a physical signal or is a noise of the apparatus that has not been properly 
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Frequency [Hz] 

FIG. 1. The Virgo sensitivity curve (from ref. [16]). 



accounted for. This might not be a great problem if the SNR is very large, but certainly with a single detector we 
cannot make a reliable detection at SNR of order one, so that the above estimates (which have been obtained setting 
SNR=1) are really overestimates. 

Virgo- Virgo. We now consider the sensitivity that could be obtained at Virgo if the planned interferometer were 
correlated with a second identical interferometer located at a few tens of kilometers from the first, and with the same 
orientation^] This distance would be optimal from the point of view of the stochastic background, since it should be 
sufficient to decorrelate local noises like, e.g., the seismic noise and local electromagnetic disturbances, but still the 
two interferometers would be close enough so that the overlap function does not cut off the high frequency range, as 
it happens, instead, correlating the two LIGOs. 

Let us first give a rough estimate of the sensitivity using h c (f), h n (f). From fig. 1 we see that we can take, for our 
estimate, hf ~ 10~ 22 Hz -1 / 2 over a bandwidth A/ ~ 1 kHz. Using T = 1 yr, eq. ( p3| ) gives 

a »«-«^»- 21 (toS5) 1/J (iF^)' < 28 > 

Requiring for instance SNR=1.65 (this corresponds to 90% confidence level; a more precise discussion of the statistical 
significance, including the effect of the false alarm rate can be found in ref. [0) gives an estimate for the minimum 
detectable value of /ioSl gw (/), 

« B W~ 3xl0 "'(ii4) a (lF4£^) ' (2S) 

This suggests that correlating two Virgo interferometers we can detect a relic spectrum with /iQJ7 gw (100Hz) ~ 3 x 10~ 7 
at SNR=1.65, or 1 x 10" 7 at SNR=1. Compared to the case of a single interferometer with SNR=1, eq. (p7|), we 



^Correlations between two interferometers have already been carried out using prototypes operated by the groups in Glasgow 
and at the Max Planck Institute for Quantum Optics, with an effective coicident observing period of 62 hours ^o|. Although 
the sensitivity of course is not yet significant, they demonstrate the possibility of making long-term coincident observations 
with interferometers. 
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gain five orders of magnitude. As already discussed, to obtain a precise numerical value one must however resort to 
eq. (0) . This involves an integral over all frequencies, (that replaces the somewhat arbitrary choice of A/ made above) 
and depends on the functional form of /igfi gw (/). If for instance ^o^gw(/) is independent of the frequency, using the 
numerical values of hf plotted in fig. 1 (see ref. fl6|| ) and performing the numerical integral we get for the minimum 
detectable h%p, gw (we give the result for a generic value of the SNR and of the integration time) 

/ SNR^ 2 1/2 



htfl™ ~ 2 x 10- 7 I— j ^-jr) » gw (/)= const.). (30) 

We see that this number is quite consistent with the approximate estimate (^), and with the value 2 x 10~ 7 reported 
in ref. pl[ . Stretching the parameters to SNR=1 (68% c.l.) and T — 4 years, the value goes down at (3 — 4) x 10~ 8 . It 
is interesting to note that the main contribution to the integral comes from the region / < 100 Hz. In fact, neglecting 
the contribution to the integral of the region / > 100 Hz, the result for /igfi™ 11 changes only by approximately 2%. 
Also, the lower part of the accessible frequency range is not crucial. Restricting for instance to the region 20 Hz 
< f < 200 Hz, the sensitivity on /igilgw degrades by less than 1%, while restricting to the region 30 Hz < / < 100 
Hz, the sensitivity on h%fil gw degrades by approximately 10%. Then, from fig. 1 we conclude that by far the most 
important source of noise for the measurement of a flat stochastic background is the thermal noise. In particular, the 
sensitivity to a stochastic background is limited basically by the mirror thermal noise, which dominates in the region 
40 Hz</<200 Hz, while the pendulum thermal noise dominates below approximately 40 Hz.p*] 

The sensitivity depends however on the functional form of fl gw (f). Suppose for instance that in the Virgo frequency 
band we can approximate the signal as 

figw(/) = ttgw(lkHz) (j^j . (31) 

For a = 1 we find that the spectrum is detectable at SNR=1.65 if ft.§f2 gw (lkHz) ~ 3.6 x 10~ 6 . For a = — 1 we 
find (taking / = 5Hz as lower limit in the integration) ft,Q^gw(lkHz) ~ 6 x 10~ 9 . Note however that in this case, 
since a < 0, the spectrum is peaked at low frequencies, and /igf2 gw (5Hz) ~ 1 x 10~ 6 . So, both for increasing or 
decreasing spectra, to be detectable ft-p^gw must have a peak value, within the Virgo band, of order a few xl0~ 6 in 
the case a — ±1, while a constant spectrum can be detected at the level 2 x 10 -7 . Clearly, for detecting increasing 
(decreasing) spectra, the upper (lower) part of the frequency band becomes more important, and this is the reason 
why the sensitivity degrades compared to flat spectra, since for increasing or decreasing spectra the maximum of the 
signal is at the edges of the accessible frequency band, where the interferometer sensitivity is worse. 

LIGO-LIGO. The two LIGO detectors are under construction at a large distance from each other, ~ 3000 km. 
This choice optimizes the possibility of detecting the direction of arrival of GWs from astrophysical sources, but it is 
not optimal from the point of view of the stochastic background, since the overlap functions cuts off the integrand in 
eq. (p) around 60 Hz. The sensitivity to a stochastic background of the LIGO-LIGO detectors has been computed in 
refs. The result, for fi gw (/) independent of /, is 

h 2 n gw ~ (5 - 6) x 10~ 6 (32) 

for the initial LIGO. The advanced LIGO aims at 5 x 10 -11 . These numbers are given at 90% c.l. in Q, and a 
detailed analysis of the statistical significance is given in Q . 

Resonant masses and Resonant mass-Interferometer. Resonant mass detectors includes bars like NAU- 
TILUS, EXPLORER and AURIGA (see e.g. refs. for reviews). Spherical HH] and truncated icosahedron 
(TIGA) pjl resonant masses are also being developed or studied. The correlation between two resonant bars and 



"Note also that it is not very meaningful to give more decimal figures in the minimum detectable value of h^Q™ 11 . Apart from 
the various uncertainties which enter the computation of the sensitivity curve shown in fig. 1, a trivial source of uncertainty is the 
fact that the computation of the thermal noises are performed using a temperature of 300K. A 5% variation, corresponding to an 
equally plausible value of the temperature, gives a 5% difference in hg£l™ n . Quoting more figures is especially meaningless when 
the minimum detectable /lo^gw is estimated using the approximate quantities h c (f),h n (f), i.e. approximating the integrand 
of eq. (Q) with a constant over a bandwidth A/. For a broadband detector these estimates typically give results which agree 
with the exact numerical integration of eq. (ph at best within a factor of two. 
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between a bar and an interferometer has been considered in refs. E6 21,2q-ES|. The values quoted in ref. [M are 



as follows (using 1 year of integration time and SNR=1; the mimimum detectable value of /igfigw grows as (SNR) , 
so the minimum detectable value at SNR=1.65 is about a factor of 3 larger). Correlating the AURIGA and NAU- 
TILUS detectors, with the present orientation, we can detect /igf2 gw — 2 x 10~ 4 . Reorienting the detectors for 
optimal correlation, which is technically feasible, we can reach tiffin ~ 8 x 1CU 5 . For AURIGA- VIRGO, with present 
orientation, /igf2 gw ~ 2 x 10~ 4 , and for NAUTILUS- VIRGO, /igf2 gw ~ 3.5 x 1CT 4 . A three detectors correlation 
AURIGA- NAUTILUS- VIRGO, with present orientation, would reach /igf2 gw ~ 1 x 10~ 4 , and with optimal orienta- 
tion /igi7 gw ~ 6 x 10 -5 . Although the improvement in sensitivity in a bar-bar-interferometer correlation is not large 
compared to a bar-bar or bar-interferometer correlation, a three detectors correlation would be important in ruling 
out spurious effects p0| . Preliminary results on a NAUTILUS-EXPLORER correlation, using 12 hours of data, have 
been reported in pffi and give a bound /igf2 gw ~ 120. 

Using resonant optical techniques, it is possible to improve the sensitivity of interferometers at special values of the 
frequency, at the expense of their broad-band sensitivity. Since bars have a narrow band anyway, narrow-banding the 
interferometer improves the sensitivity of a bar-interferometer correlation by about one order of magnitude |2lf| . 

While resonant bars have been taking data for years, spherical detectors are at the moment still at the stage of 
theoretical studies (although prototypes might be built in the near future), but could reach extremely interesting 
sensitivities. In particular, two spheres with a diameter of 3 meters, made of A15056, and located at the same site, 
could reach a sensitivity /igfi gw ~ 4 x 10~ 7 fTofl . This figure improves using a more dense material or increasing the 
sphere diameter, but it might be difficult to build a heavier sphere. Another very promising possibility is given by 
hollow spheres The theoretical studies of ref. p3 suggest a minimum detectable value /igfi gw ~ 10~ 9 at / = 200 
Hz. 

IV. THE ENERGY SCALES PROBED BY RELIC GRAVITATIONAL WAVES 

Let us consider the standard Friedmann- Robertson- Walker (FRW) cosmological model, consisting of a radiation- 
dominated (RD) phase followed by the present matter-dominated (MD) phase, and let us call a(t) the FRW scale 
factor. The RD phase goes backward in time until some new regime sets in. This could be an inflationary epoch, 
e.g. at the grand unification scale, or the RD phase could go back in time until Planckian energies are reached and 
quantum gravity sets in, i.e., until t ~ tp\ ~ 5 x 10~ 44 s. If the correct theory of quantum gravity is provided 
by string theory, the characteristic mass scale is the string mass which is somewhat smaller than the Planck mass 
and is presumably in the 10 17 — 10 18 GeV region, and the corresponding characteristic time is therefore one or two 
orders of magnitude larger than tpi. The transition between the RD and MD phases takes place at t = t eq , when the 
temperature of the Universe is of the order of only a few eV, so we are interested in graviton production which takes 
place well within the RD phase, or possibly at Planckian energies. 

A graviton produced with a frequency /* at a time t = t*, within the RD phase has today (t — to) a red-shifted 
frequency /o given by fo = /*a(i*)/a(to)- To compute the ratio a(t*)/a(to) one uses the fact that during the 
standard RD and MD phases the Universe expands adiabatically. The entropy per unit comoving volume is S = 
const. gs(T)a 3 (t)T 3 , where gs{T) counts the effective number of species §]. In the standard model, at T>300 GeV, 
g s {T) becomes constant and has the value g s (T*) = 106.75, while today g S (T ) ~ 3.91 f| and T = 2.728±0.002K |§. 
Using g S (T*)a 3 (U)T 3 = g s (T )a 3 (t Q )T 3 one finds Q 

/0 = / ,f), 8 .0x 1 0->V.( "° (33) 
a(*o) \gs{T*)J \ T* J 

The first point to be addressed is what is the characteristic value of the frequency /* produced at time f » , when the 
temperature was T*. One of the relevant parameters in this estimate is certainly the Hubble parameter at time of 
production, H(t*) = H*. This comes from the fact that H^ 1 is the size of the horizon at time t*. The horizon size, 
physically, is the length scale beyond which causal microphysics cannot operate (see e.g. ||, ch. 8.4), and therefore, 
for causality reasons, we expect that the characteristic wavelength of gravitons or any other particles produced at 

smaller.^] 



^On a more technical side, the deeper reason has really to do with the invariance of general relativity under coordinate 
transformations, combined with the expansion over a fixed, non-uniform, background. Consider for instance a scalar field <f>(x) 
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Therefore, we write A* = eH~ x . The above argument suggests e < 1. During RD, H% = (8n/3)Gp la d- The 
contribution to p from a single species of relativistic particle with gi internal states (helicity, color, etc.) is p ra d = 
gi(w 2 / 3Q)T 4 for a boson an p = (7/8).gi(7r 2 /30)T 4 for a fermion. Taking into account that the i-th species has 
in general a temperature Tj 7^ T if it already dropped out of equilibrium, we can define a function g(T) from 
Prad = (tt 2 /3(%(T)T 4 . Then § 

9(T)= E ^(|) 4 + ^ E (34) 

i— bosons z— fcrmions 

The sum runs over relativistic species. This holds if a species is in thermal equilibrium at the temperature TV If 
instead it does not have a thermal spectrum (which in general is the case for gravitons) we can still use the above 
equation, where for this species Ti does not represent a temperature but is defined (for bosons) from pi = g i (7r 2 /30)T l 4 , 
where pi is the energy density of this species. The quantity gs(T) used before for the entropy is given by the same 
expression as g(T), with (T^/T) 4 replaced by (Ti/T) 3 . We see that both g(T) and gs(T) give a measure of the effective 
number of species. For most of the early history of the Universe, g{T) = gs(T), and in the standard model at T> 
300 GeV they have the common value g* = 106.75, while today g(To) — 3.36, gs{To) = 3.91 Therefore 



H * ~ 90M 2 , ' [6b) 



and, using /» = i?*/e, eq. ( p3| ) can be written as 31 



7 1 (%\ /^V- 



/0 - L65Xl0W UlGivJ^» (36) 

This simple equation allows to understand a number of important points concerning the energy scales that can be 
probed in GW experiments. The simplest estimate of /* corresponds to taking e = 1 in cq. ( |36| ) [TT] . In this case, 
we would find that a graviton observed today at the frequency fo = 1Hz was produced when the Universe had a 
temperature T t ~ 6 x 10 6 GeV. Using the relation between time and temperature in the RD phase, 



2.42 ' A ' 1 ' " " 



1/2 

g* 



sec , (37) 



we find that this corresponds to a production time t* ~ 7 x 10 21 sec, and at time of production this graviton had 
an energy E, ~ 0.3 MeV. At f = 100Hz we get t* ~ 7 X 10~ 25 sec, T t ~ 6x 10 s GeV and = 3 GeV. These would 
be, therefore, the scales relevant to Virgo and LIGO. For a frequency fo = 10 -4 Hz, relevant to LISA, we would get 
instead U=7 x 10" 13 sec, T* - 600 GeV. 

However, the estimate A* ~ H^ 1 , or e ~ 1, can sometimes be incorrect even as an order of magnitude estimate. In 
the Appendix we illustrate this point with two specific examples, one in which the assumption A* ~ H^ 1 turns out 
to be basically correct, and one in which it can miss by several orders of magnitudes. Both examples will in general 
illustrate the fact that the argument does not involve only kinematics, but also the dynamics of the production 
mechanism. 

From eq. ( |36"1) we see that the temperatures of the early Universe explored detecting today relic GWs at a frequency 
fo are smaller by a factor approximately equal (for constant t/») to e, compared to the estimate with e = 1. Equivalently, 



and expand it around a given classical configuration, <j>{x) = 4>o(x) + S(j>(x). Under a general coordinate transformation x — > x', 
by definition a scalar field transforms as <j>{x) — * <f>'(x') — 4>(x). However, when we expand around a given background, we keep 
its functional form fixed and therefore under x —* x' , <f>o(x) — + <f>o(x'), which for a non-constant field configuration, is different 
from <f>o(x)- It follows that the perturbation 5(j>{x) is not a scalar under general coordinate transformations, even if (f>(x) was a 
scalar. The effect becomes important for the Fourier components of 8cj>{x) with a wavelength comparable or greater than the 
variation scale of the background 4>q(x). (We are discussing a scalar field for notational simplicity, but of course the same holds 
for the metric tensor g^v). In a homogeneous FRW background the only variation is temporal, and its timescale is given by the 
H~ . Therefore modes with wavelength greater than H~ l are in general plagued by gauge artefacts. This problem manifests 
itself, for instance, when computing density fluctuations in the early Universe. In this case one finds spurious modes which can 
be removed with an appropriate gauge choice, see e.g. ref. EL sect. 9.3.6 or ref. |32j. 
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a signal produced at a given temperature could in principle show up today in the Virgo/LIGO frequency band 
when a naive estimate with e — 1 suggests that it falls at lower frequencies. 

There is however another effect, which instead gives hopes of exploring the Universe at much higher temperatures 
than naively expected, using GW experiments. In fact, the characteristic frequency that we have discussed is the value 
of the cutoff frequency in the graviton spectrum. Above this frequency, the spectrum decreases exponentially |53| , and 
no signal can be detected. Below this frequency, however, the form of the spectrum is not fixed by general arguments. 
Thermal spectra have a low frequency behaviour Q gw (f) ~ f 3 , as we read from eq. ( pl| ) inserting a Bose-Einstein 
distribution for rif, so that, at low /, rif ~ 1//. However, below the Planck scale gravitons interact too weakly to 
thcrmalizc, and there is no a priori reason for a ~ / 3 dependence. The gravitons will retain the form of the spectrum 
that they had at time of production, and this is a very model dependent feature. However, from a number of explicit 
examples and general arguments that we will discuss below, we learn that spectra flat or almost flat over a large range 
of frequencies seem to be not at all unusual. 

This fact has potentially important consequencies. It means that, even if a spectrum of gravitons produced during 
the Planck era has a cutoff at frequencies much larger than the Virgo/LIGO frequency band, still we can hope to 
observe in the WHz-lkHz region the low- frequency part of these spectra. In the next subsection we will therefore 
discuss what signals can be expected from the Universe at extremely high (Planckian, string or GUT) temperatures. 



V. TOWARD THE PLANCK ERA? 



The scale of quantum gravity is given by the Planck mass, related to Newton constant by G — l/M™. More 
precisely, since in the gravitational action enters the combination 8irG, we expect that the relevant scale is the 
reduced Planck mass M = M p1 /(8tt) 1 / 2 ~ 2.44 x 10 18 GeV. Using eq. @ with T* = M and e = 1 gives 

fo ~ 400 (^) V6 GHz , (T, - Mn/Vto) . (38) 

The dependence on is rather weak because of the power 1/6 in eq. J36|) . For g* = 1000, fo increases by a factor 
~ 1.5 relative to g* = 100. For T* = A/qut ~ 10 16 GeV, and ~ 220 (likely values for a supersymmetric unification), 

fo ~ 2 GHz , (T, = M GUT ) • (39) 

Using instead the typical scale of string theory, Mg, the characteristic frequency is between these two values, since 
Mg is expected to be approximately in the range 10 _2 <M5/Mpi<10 _1 |34|| . Actually, in the estimate of fo, precise 
numbers depend on details of the criteria used. For instance, in the context of string cosmology [ ^5|j36| it might be 
more meaningful to perform the estimate setting iJ» ~ Ms rather than T* ~ Ms, since the Hubble constant during the 



large curvature regime is fixed by stringy corrections, whose scale is given by Ms |p7|| . In this case, from eq. (35), the 
estimate for T* is T* - (90/87r 3 5») 1 / 4 (Mp 1 Ms) 1/2 instead of T* ~ M s . Thus, the dependence of f on 5» is different 
(~ <?* X ^ 12 instead of the factor gl^ 6 that appears in eq. (^8|)), and there appears also a dependence on M$/Mp\ [fS6"[ . 
However, numerically the difference is not very significant, since if say g* — 100, then (QO/Stt 3 ^) 1 / 4 ~ 0.25 while for 
<7* = 1000, (QO/Stt 3 ^) 1 / 4 ~ 0.14, and this number is partially compensated by (Mpj/M^) 1 / 2 , that ranges between 
approximately 10 and 3; thus, for fo the overall number with the estimate T* ~ Ms or with ~ Ms can differ, say, 
by a factor of two, which is anyway beyond the accuracy of these order of magnitude estimates. 

The typical frequency region for GUT/string/Planck physics is therefore between the GHz and a few hundreds 
GHz. This is very far from the region accessible to interferometers, which arrive at most to 10 kHz, see fig. 1. 

One might consider the possibility of building a detector in the GHz region to search for a signal from the Planck 
era. GW detectors based on microwave superconducting cavities that could operate at these frequencies have indeed 
been studied in the literature p8| , p9[ in the years 1978-1980. However, from eq. @ we see that a level /ip^gw ~ 10 6 
corresponds, at / = 1 GHz, to h c ~ 10~ 30 , which is many orders of magnitude below the sensitivities discussed in Q. 
Still, given the great interest that a detector in this frequency range could have, it could be worthwhile to investigate 
them further Q 



'"''These cavities could also be constructed so that they operate in the MHz region [^0|, where the requirements on h c are less 
stringent, and the wavelength is long enough to perform a correlation between two detectors. A relic background detected at 
these frequencies would be extremely interesting since it would be unambiguously of cosmological origin, see sect. 7B. 
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It is clear from the above estimates of the characteristic frequency, that a necessary condition for observing at Virgo 
or LIGO a signal from the Planckian era is that we have a relic gravitational spectrum which is not peaked around 
the cutoff frequency, but rather is almost flat from the GHz region, where it has the cutoff, down to at least the kHz 
region, where interferometers can operate. 

A thermal spectrum behaves as /igf2g W ~ / 3 and the value in the kHz region would be a factor ~ 10 -18 smaller than 
a peak value at the GHz, and therefore would be obviously unobservable. However, as already remarked, gravitons do 
not thermalize and there is no a priori reason for a ~ / 3 spectrum. On the contrary, there are good general reasons 
and various explicit examples favoring almost flat spectra, at least over some range of frequencies. 

A simple and instructive example has been given by Krauss El]. Consider a global phase transition in the early 
Universe, associated with some scalar field <p that, below a critical temperature, gets a vacuum expectation value 
((j)). For causality reasons, in an expanding Universe the scalar field cannot have a correlation length larger than the 
horizon size. So, even if the configuration which minimizes the energy is a constant field, the field will be constant 
only over a horizon distance. During the RD phase the horizon expands, and the field will relax to a spatially uniform 
configuration within the new horizon distance. This relaxation process will in general produce gravitational waves, 
since there is no reason to expect (T^ v ) for the classical field just entering the horizon to be spherically symmetric. 
A simple estimate of the spectrum produced makes use of the fact that the characteristic scale of spatial or temporal 
variation of the field is given by the inverse of the Hubble constant, H~ x , at the value of time under consideration. 
Then, the energy density of the field is p ~ (d<j)) 2 ~ (<f)) 2 H 2 . The total energy in a Hubble volume is pH~ 3 ~ (0) 2 iJ _1 
and the corresponding quadrupole moment is Q ~ (pH~ 3 )H~ 2 ~ (cf>) 2 H~ 3 , times a constant smaller than one, which 
measures the non-sphericity of the configuration, and which is not very relevant for order of magnitude estimates of 
the frequency dependence. The energy liberated in GWs in a horizon time is then given by fill 

AE ~ H- 1 x Luminosity ~ H^G (^jjjf) ~ G(0) 4 iJ" 1 , (40) 



and the energy density of GW produced is 

dp gw AE 
dlogf ~ IF 3 



G{cj)) 4 H 2 . (41) 



This is the energy density at a frequency / ~ H, at a value of time t = t(f) which is the time when the mode with 
frequency / enters the horizon. So, dp/ dlogf in the above equation is evaluated at different values of time for different 
frequencies. However, in the RD phase, the energy density scales like p ~ l/a 4 (t) ~ t~ 2 and H ~ i" 1 so that p scales 
as H 2 . Therefore, eq. ( ftlj ) means that, if we evaluate dp/dlogf at the same value of time, for all frequencies which 
at this value of time are inside the horizon, the spectrum is flat, fto^gw(/) ~ const. This is a rather nice example of 
how a flat spectrum can follow simply from the requirement of causality and simple dimensional estimates. 

The best known example of an (almost) flat spectrum is the amplification of quantum fluctuations in the inflation- 
RD phase transition [^| 16| . In this case, during the inflationary phase the horizon size H~ x is approximately constant 
and perturbations exit from the horizon, since the physical wavelength scales as ~ a(t). The spectrum of quantum 
fluctuations for gravitational waves during inflation is (see e.g. M, pg.289) Ah + _ x (f) ~ H/Mj>\, where H is the 
Hubble constant at time of crossing. Different wavelengths cross at different time, but since H does not vary during 
inflation, this is really independent of /. When the mode is outside the horizon, its amplitude stays constant, again 
for causality reasons. During the RD phase the horizon size grows faster than a(t) and so finally the perturbation 
reenters the horizon, and its energy density appears in the form of gravitons. Its energy density as it crosses back is 
dp/dlogf ~ h 2 r~j f 2 {Ah) 2 . As in the previous example, this is the value of dp/dlogf at frequency /, at the value of 
time when this wavelength crosses back inside the horizon, so it is evaluated at different times for different frequencies, 
and a dependence ~ f 2 means that the spectrum, when dp/dlogf is evaluated at the same time for all frequencies, 
is flat. So, again in this example a basic ingredient for a flat spectrum was causality, togheter with the fact that H is 
constant during inflation. 

The same mechanism of amplification of vacuum fluctuations produces an interesting stochastic background in 
string cosmology. In recent years Gasperini and Veneziano have proposed an inflationary model based on string 
theory ]35|]. In this model the Universe starts at low curvatures; at this stage it is described by the lowest order 
effective action of string theory. Then it evolves toward a regime of large curvature, where higher order corrections to 
the string effective action stabilize the growth of the curvature p7|j4^ ], which is now of order one in string units, and 
the evolution should then be matched to the standard RD phase |48f! The big-bang, i.e. a state of the Universe with 
Planckian energies and curvature, rather than being a postulated initial condition, is the end result of an evolution 
(the 'pre-big-bang' phase) which started from a cold Universe at low curvature. 
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The form of the spectrum of relic GWs has been discussed in various recent papers |49|- |54| , p6l ancl depends on the 
rate of growth of the curvature and of the dilaton field in the various phases of the model. Very low frequencies are 
sensitive to the behavior in the low curvature pre-big-bang regime. In this regime both the curvature and the dilaton 
are growing and the spectrum turns out to grow as /igf2 gw ~ / 3 log 2 /. In the large curvature regime, according to a 
possible scenario p7j] the Hubble parameter is constant, while the dilaton can grow, with a derivative parametrized 
by an unknown constant. In the limit in which the derivative of the dilaton in this phase goes to zero, the spectrum 
in the region f s < f < / pca k is approximately flat. The value of f s is determined by the duration of this intermediate 
large curvature phase (/ pca k//s is in fact the total red-shift during this phase), and it is not calculable without a 
detailed knowledge of stringy phenomena. With a sufficiently long string phase, it is possible to have f s smaller than, 
say, 1kHz. The cutoff value / pca k is instead in the region of a few hundred GHz, as already discussed. Therefore, 
if in the large curvature regime we have a long inflationary phase with an almost constant dilaton, we have again a 
spectrum which is practically flat over a very large range of frequencies.^ 

Another example of relic GW spectrum which is almost flat over a huge range of frequencies is provided by GWs 
produced by the decay of cosmic strings Js^j . In this case the spectrum has a peak around / ~ 1CP 12 Hz, where 
/igilgw can reach a few times 10 -6 , and then it is almost flat over a huge range of frequencies, from / ~ 10 -8 Hz to 
the GHz region, where it has the cutoff (fixed by the arguments previously discussed, considering that the relevant 
scale for cosmic string is MIqut)- The basic reason for such a behavior is the scaling property of the string network, 
which says that a single lengthscale, the Hubble length, characterizes all properties of the string network (5^ Jll[] . The 
network of strings evolves toward a self-similar configuration, with small loops being chooped off very long strings, 
and the typical radiation emitted by a single loop has a wavelength related to the length of the loop. 

In conclusion, from these examples we learn that GW spectra almost flat over a large range of frequencies are not 
at all unusual in cosmology, and often the approximate flatness of the spectrum is a consequence of rather general 
principles, like causality, or the existence of a single lenghtscale in the problem. Therefore, there is some hope that 
we need not go in the GHz region to explore Planck scale physics. 



VI. THE EFFECTS OF INFLATION 



In spite of the fact that the standard cosmological model based on a RD and a MD phase is not contradicted by 
any experimental data, inflation (i.e., a stage of exponential expansion, or more generally a period characterized by 
the condition a(t) > on the FRW scale factor a(t)) has many compelling features, see e.g. refs. Q,|5(||37| for reviews. 
At present, there is no 'standard model' of inflation. The mechanism is very general, and can be implemented in 
different models. In the models which have been 'traditionally' studied, inflation is driven by the potential energy of 
a scalar field. 

Independently of the details, these models predict an almost flat gravitational wave spectrum extending from some 
high-frequency cutoff down to extremely small frequencies, / ~ 10~ 16 Hz, and rising as l// 2 as we go at even lower 
frequencies, such as / ~ 10 -18 Hz, that corresponds to wavelengths of the order of the present Hubble distance 
i? _1 (to). The value of ft,Qf2 gw (/) for this spectrum depends on the Hubble constant during inflation, i/r>S) which in 
turn is related to the energy scale Mj n fl at which inflation takes place, i?DS ~ Miaa/Mpi. A strong signal in GWs at 
wavelengths A ~ if _1 (£o) would produce a stochastic red-shift on the frequencies of the photons of the 2.7K radiation, 
and a corresponding fluctuation 8T/T in their temperature. The measurements of ST/T from COBE therefore put 
an upper bound on the intensity of the relic GWs produced at these long wavelengths, that translates into an upper 
bound on the Hubble constant during inflation and on the scale Mj n fl where inflation occurs |4^j5q| . Using COBE 
four-year data (and taking properly into account the relation between the amplitude at low frequency and the tilt of 
the spectrum, which is not exactly flat) gives an upper bound |p8| , p7| 

M infl < 3.4 x 10 16 GeV . (42) 

A crucial consequence of this bound would be that we cannot see today a cosmological signal from the Planck scale. 
Every particle production which took place at Planckian times has been exponentially diluted by the subsequent 



§ § Note that, while in standard inflationary models a long inflationary phase at large curvature would provide too much GWs 
at the frequecies relevant for the COBE bound (see below), this is not the case for string cosmology because of the behaviour 
~ f 3 log 2 / that sets in at low /. 
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inflationary stage. In this case the maximum energy that can be explored by cosmological observations is the reheating 
temperature after inflation, which is equal to Af; n fi in the case of a perfectly efficient reheating, and is in general lower. 

However, the above conclusion is not unescapable. We have already considered above the inflationary model based 
on string theory proposed by Gasperini and Veneziano ]35|] . In this model we have automatically an inflationary 
evolution during the pre-big-bang phase, triggered by the kinetic energy of the dilaton field. The bound ( |42| ) in this 
case does not apply because the spectrum of relic gravitational waves is not flat over the whole frequency range, from 
COBE frequencies / ~ 1CP 16 up to the cutoff frequency, as in standard inflation, but rather goes like ~ / 3 log 2 / at 
sufficiently small /, thereby automatically avoiding the COBE bound by many orders of magnitude G9[. 



Refs. [p5 59 1 discuss whether this inflationary phase can solve the kinematical problems of standard cosmology 



without fine-tuned initial conditions. In this paper we do not want to enter a discussion on the positive and negative 
aspects of any specific model. We rather take the model of refs. |55| as an indication that the bound ( [l2| ) is not 
necessarily unescapable, and we discuss in subsections A and B two different scenarios which depend on whether 
this bound is obeyed or not. The amplification of vacuum fluctuations is common to both scenarios and is discussed 
in subsection C, where we also consider the possibility of combining the pre-big-bang scenario with potential-driven 
inflation. 



A. Cosmological signals from the reheating era 

If we do have an inflationary phase at a scale Mj n fl that satisfies the bound (^), there is no gravitational analog 
of the 2.7K radiation. At the Planck scale photons and gravitons can be produced simply by thermal collisions, with 
probably similar production rates. However, the particles produced in this way at the Planck (or string) scale are 
exponentially diluted by the subsequent inflationary phase, and the photons that we see today in the cosmic microwave 
background radiation (CMBR) have been produced during the reheating era which terminated the inflationary phase. 
The highest energy scale that can be probed by cosmological observations in this case is given by the reheating 
temperature (the amplification of vacuum fluctuations can provide an exception to this statement, see subsect. C). 

Since the reheating temperature T r h < M- m & <C Mpi, thermal collisions are by now unable to produce a substantial 
amount of gravitons. However, in this case relic GWs can be produced through some non-equilibrium phenomenon 
connected with the reheating process. For a process that takes place at the reheating temperature T r h the estimate of 
the characteristic frequency is given by eq. ([36j) with T* = T r h. In principle the reheating temperature can be between 
a minimum value at the TeV scale (since this is the last chance for baryogenesis, via the anomaly in the electroweak 
theory) and a maximum value T r h ~ Mj n fl<3 x 10 16 GeV. In supersymmetric theories, the gravitino problem give a 
further constraint T r h^l0 9 GeV p|]6C|]. More precisely, ref. (6^] gives a bound on T r h ~ 10 6 — 10 9 GeV for a gravitino 
mass 100 GeV-1 TeV, and a bound 10 11 — 10 12 GeV independent of the gravitino mass. 

So, phenomena occurring at reheating will manifest themselves with cutoff frequencies between 10 -4 and 10 9 Hz, 
if we set e = 1 in eq. (|36|). In typical non-supersymmetric models the reheating temperature is large, say 10 14 GeV, 
corresponding to /q ~ (l/e)10 7 Hz, and in this case interferometers or resonant masses could only look for their 
low-frequency tails. If instead the rehating temperature is T r h ~ 10 9 GeV, non-equilibrium phenomena taking place 
during reheating would give a signal just in the LIGO/Virgo frequency band. 

A promising mechanism for GW generation during reheating is bubbles collision, in the case when inflation ter- 
minates with a first order phase transition |]6l| , ^T| . In this case the relevant estimate is /o ~ (l/e)10 7 Hz. One 
should however keep in mind the possibility that bubble nucleation occurs before the end of inflation, and then the 
cutoff frequency would be red-shifted by the subsequent inflationary evolution toward lower values, possibly within 
the Virgo/LIGO frequency range p^fl . 

The other possibility is that reheating occurs through the decay of the inflaton field. In this case, in a very general 
class of models, there is first an explosive stage called preheating, where the inflaton field decays through a non- 
perturbative process known as parametric resonance J63| , and at this stage there are mechanisms that can produce 
GWs, see ^,^|. We believe that in these cases a reliable estimate of the characteristic frequency is really very difficult 
to obtain, since the relevant parameter, which is the value of the Hubble constant at time of production, depends on 
the complicated dynamics of preheating and on the specific inflationary model considered. Ref. ]64|] examines models 
where the characteristic frequency turns out to be in the region between a few tens and a few hundreds kHz. 



15 



B. Cosmological signals from the Planck era 



The second possibility is that there is no inflationary phase after the Planck era. This could happen if there is no 
inflationary phase at all (although in this way we would lose the advantages of inflation, so that this possibility does 
not seem theoretically very appealing), or if inflation is realized in the 'pre-big-bang' phase of string cosmology, as 
suggested in [ |35|| . Independently of other considerations, the pre-big-bang model is interesting in this context because 
it is an explicit example of the fact that physics at the Planck or string scale can invalidate the considerations leading 
to the bound ([f2|), without necessarily giving up the idea of inflation. 

If there is no inflationary phase after the Planck era, the photons of the CMBR have been produced at the Planck 
(or string) scale, when gravitons are also produced with probably similar rates, and therefore we expect a GW 
analog of the 2.7 K radiation. At Planckian times both gravitons and photon will have approximately the same 
characteristic energy. Then, gravitons decouple and their characteristic energy scales like l/a(t), while photons are 
in thermal equilibrium and their temperature evolves so that the entropy ~ gs(T)T 3 a 3 (t) is constant. Therefore the 
characteristic energy of the gravitons today, Tq, is related to the present temperature of the photon, To = 2.728±0.002, 

by | 

3.9l\ 1/3 m „„„/100^ 1/3 



T G = — T ~ 0.93 K , (43) 

V 9* ) \ 9* J 

where we used g s (T ) ~ 3.91 §] and is the effective number of degrees of freedom at time of production. As we 
already discussed, in the case of the gravitons Tq should not be considered a temperature, since there is no a priori 
reason why gravitons should have a thermal spectrum, but rather is a characteristic energy, close to the high-energy 
cutoff in the spectrum. The corresponding frequency is /o — 120 (100/g*) 1 / 3 GHz; this is the cutoff of a spectrum 
which extends toward lower frequencies. If its behaviour is, say, /ig^gw ~ f 3 there is no hope of observing it in the 
kHz region, while if it is flatter it can be more interesting for present experiments. 

Actually, since in this case we are looking for a signal coming from the quantum gravity regime, the situation can 
be richer. For instance, if the correct theory of quantum gravity is provided by string theory, it may happen that 
massive string modes are highly excited in the large curvature regime, as indeed happens in the scenario discussed in 
ref. (4?]]. In this case the decay of excited string modes can be an interesting mechanism for production of GWs. In 
fact, for large level number N the decay width of the string is dominated by transitions of the form N — > TV — 1 with 
emission of a graviton or another massless particle [Q. The mass levels of the string depend on N as Mjy ~ N so 
that, in string units, AM N ~ l/\/~N, and therefore these decays would produce a set of peaks at today frequencies 
In = fo/VN. For large N these frequencies are close to each other and it is possible that the width will be large, 
so that we will not have lines but rather a continuous spectrum, but in any case this mechanism shifts some of the 
energy from the GHz region toward a more accessible frequency range. 

Note also that the value f ~ 120(100/5*) 1/3 GHz is really only an upper bound, since the production through 
decay of excited string levels would take place before the onset of the RD phase, and a further red-shift should be 
taken into account. This cannot be computed in the absence of a detailed knowledge of the cosmological evolution in 
the large curvature regime, but in general it might be possible that the effect of the redshift and/or the factor 1/yN 
shift this radiation in an interesting frequency band M7J. 



C. Amplification of vacuum fluctuations 

A very general and well studied mechanism for the production of relic GWs is the amplification of vacuum fluctua- 
tions due to a non-adiabatic transition between different regimes for the evolution of the scale factor. The mechanism 
is common to the scenarios discussed in the two previous subsections. 

A well known example is the transition between a (potential-driven) inflationary phase and the standard RD 
phase [p2|— p6[ . At first sight, it is not promising for detection, since once the bound (^2|) is obeyed, at LIGO/Virgo 
frequencies one gets h^Q,^ ~ 10~ 15 j67|, very far from the sensitivity of first and even of the planned second generation 
experiments. This is due to the fact that in this model /ipilgw is almost flat at all frequencies, from COBE frequencies 
up to the cutoff of the spectrum, and therefore the bound imposed by COBE at very low frequencies gives a bound 
at all higher frequencies. One should however keep in mind that these computations are performed assuming that 
before the transition the Universe is in the DeSitter invariant vacuum. If instead before the transition we have njj 
gravitons per cell of the phase space, the total number of gravitons after the transition is 
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n k =nl + 2\fo\ 2 (iil+ l ^j , (44) 

where (3k is the relevant Bogoliubov coefficient (we are considering the case fikk' — PkSkk', as is usually the case in 
cosmological backgrounds). If n? = we have amplification of the 'half quantum' of vacuum fluctuations, with an 
amplification factor 2|/3/ c | 2 . But also any preexisting non-vanishing n® is amplified. One can expect that before the 
transition <C 1 because a phase space density of order one would correspond to a large energy density, since at 
early time these modes have very high frequency f44}| . However, the argument need not be true if there has been some 
mechanism which pumped the energy into these modes, as may happen for instance with a previous amplification 
of vacuum fluctuations at the string scale, as in For a transition between a DeSitter phase (if=const.) and 

RD, the Bogoliubov coefficients depend on k as 0k ~ 1/fc 2 if before the transition — 0, then after the 



transition nk — \(3k\ ~ 1/fc ~ 1// and eq. (21) shows that the spectrum for /igft gw (/) is flat (actually, during 



slow- roll inflation, the Hubble parameter is not exactly constant, and this results in a small tilt in the spectrum, 
see e.g. Q). If instead for some fc we have n° ^ and \/3k\ 1, the dependence of /igSlgw on / is of the form 
/i§f2 gw (/) ~ n® + (1/2), since the factor / 4 in eq. ( pl| ) is cancelled by \(3k\ 2 ~ I// 4 - We see that the spectrum of 
vacuum fluctuations after an inflationary phase is determined by the occupation numbers before the onset of inflation. 
This is an example (actually, the only example of which we are aware) of the fact that an inflationary phase does not 
erase every information about the previous stage of the Universe. A measurement of a relic GW spectrum, and hence 
of rik, if the Universe underwent a stage of inflation, would allow us to reconstruct the state of the Universe (i.e. n°) 
before inflation. The technical reason behind this is that the quantity which is amplified is a number of particle per 
unit cell of the phase space, and the volume of a cell of the phase space, d 3 xd 3 fc/(27r) 3 , is unaffected by the expansion 
of the Universe, contrarily to a physical spatial volume d 3 x ~ a 3 (t). 

These considerations can be applied to the case in which the pre-big-bang model of Gasperini and Veneziano is 
followed, rather than by the standard RD cosmology, by a phase of inflation, driven by the potential energy of some 
inflaton field. Indeed, it was stressed already in the first papers on pre-big-bang cosmology ||^] that the pre-big-bang 
scenario is by no means to be regarded as an alternative to standard (even inflationary) cosmology, but rather it is 
a completion of the standard picture beyond the Planck era. While most efforts to date have been concentrated in 
matching the pre-big-bang phase to a RD phase, one can certainly consider the matching to, say, a stage of chaotic 
inflation. This could have some attractive features, since the inflationary stage that takes place during the pre- 
big-bang phase, independently of whether it is long enough to solve the horizon/flatness problems [^9|, can anyway 
prepare sufficiently homogeneous patches that are need to start chaotic inflation. The spectrum of relic GWs in a 
model with a pre-big-bang phase and a successive stage of potential-driven inflation would still retain at low / the 
form ft-o^gw ~ / 3 log 2 / obtained in string cosmology, and therefore the comparison of the relic GW spectrum with 
COBE anisotropies would give no bound on the scale M{ n Q of the potential-driven inflation. 

Another interesting example, considered in |39|, is the case in which, before the onset of an inflationary phase, the 
Universe was already in a state of thermal equilibrium, so that n9 is given by a Bose-Einstein distribution. This is a 
natural initial condition when the inflationary phase is triggered by a first order phase transition, since this requires 
a homogeneous thermal state as initial condition. This example works in the opposite direction compared to the case 
of string cosmology followed by potential inflation. In fact, now at small fc we have = (exp(Ek/T) — 1) _1 ~ 1/fc, 
and the spectrum is enhanced at low frequency, and totally neglegible at Virgo/LIGO frequencies. This also results 
in a stronger bound on the inflation scale compared to eq. ( fl2| ) . 

A different possibility is that the equation of state in the early Universe is modified, and this can result in growing 
spectra so that again we can obey the COBE bound at / ~ 1CP 16 Hz, but still we could have a large signal at higher 
frequencies. Ref. Q considers the case of a phase of the early Universe with scale factor expanding as a{r/) ~ 
where 77 is conformal time. This corresponds to an equation of state p = 7/? with 7 = (1 — /3)/(3 + 3/3). Inflation 
corresponds to 7 = — 1 and then to (3 = —2. At frequencies / > 10~ 16 Hz the amplification of vacuum fluctuations 
at the transition between this phase and the standard RD phase gives a spectrum with h c (f) ~ f /3+1 and therefore, 
from eq. ([ljj), ft-o^gw ~ f 2/3+4 , which of course reproduces a flat spectrum in the inflationary case, (3 = —2. The 
parameter (3 is related to the spectral index of density perturbations n by n = 2(3 + 5, so that (3 = — 2 corresponds 
to n = 1, the Harrison-Zeldovich spectrum. The value of n can in principle be derived from COBE observations. 
Two different groups, analyzing the same data set in different ways, have produced the results n — 1.2 ± 0.3 |7lJ| 
(therefore consistent with a Harrison-Zeldovich spectrum) and n = 1.84 ±0.29 |^2| . So, systematic errors presumably 
dominate the analysis, but there is the possibility that n > 1 and therefore (3 > — 2, which would give a growing GW 
spectrum /&o^gw(/) ~ f n l - However, if n is too large (like the value n — 1.84) one should also find a way to 
cutoff the spectrum before it grows too large and exceeds existing bounds like the nucleosynthesis bound, see sect. 7. 

If instead there is no phase of potential-driven inflation, while there is an inflationary pre-big-bang phase as suggested 
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in [3^ ], than one finds a relic GW spectrum that grows like / 3 log 2 / at small /, thereby avoiding the COBE bound, 
and which has the cutoff in the GHz region, as we already discussed. Whether the spectrum is approximately flat 
from the GHz down to the kHz region relevant for Virgo/LIGO depends on some unknown parameters of the large 
curvature phase, in particular the rate of growth of the dilaton and the duration of this phase. For some values of 
these parameters (almost constant dilaton and long large curvature phase) the spectrum is almost flat in a large range 
of frequencies and it might be relevant at the frequencies of present experiments ]49|- |54] , |36| ] . 



VII. CHARACTERISTIC INTENSITIES 

A. The nucleosynthesis bound 

Nucleosynthesis successfully predicts the primordial abundances of deuterium, 3 He, 4 He and 7 Li in terms of one 
cosmological parameter 77, the baryon to photon ratio. In the prediction enter also parameters of the underlying 
particle theory, which are therefore constrained in order not to spoil the agreement. In particular, the prediction is 
sensitive to the effective number of species at time of nucleosynthesis, = g(T ~ MeV). With some simplifications, 
the dependence on can be understood as follows. A crucial parameter in the computations of nucleosynthesis is 
the ratio of the number density of neutrons, n n , to the number density of protons, n p , As long as thermal equilibrium 
is mantained we have (for non-relativistic nucleons, as appropriate at T ~ MeV, when nucleosynthesis takes place) 
n n /n p = exp(— Q/T) where Q = m n — m p ~ 1.3 MeV. Equilibrium is mantained by the process pe <-> nv, with width 
Fpe^nvi as long as T pe ^> ni/ > H. When the rate drops below the Hubble constant H, the process cannot compete 
anymore with the expansion of the Universe and, apart from occasional weak processes, dominated by the decay of 
free neutrons, the ratio n n /n p remains frozen at the value exp(— Q/Tf), where Tf is the value of the temperature at 
time of freeze-out. This number therefore determines the density of neutrons available for nucleosynthesis, and since 
practically all neutrons available will eventually form 4 He, the final primordial abundance of 4 He is very sensitive to 
the freeze-out temperature Tf. Let us take for simplicity r ; „ — „„ — G F T (which is really appropriate only in the 
limit T > Q)- The Hubble constant is given by H 2 = (8ir/3)Gp, where p include all form of energy density at time 
of nucleosynthesis, and therefore also the contribution of primordial GWs. As usual, it is convenient to write the 
total energy density p in terms of see eqs. (p^), as p — (n 2 / '30) <?*T 4 . We recall that, for gravitons, the quantity Tj 
entering eq. j34| ) is defined by p gw = 2{ir 2 /30)Tf, and this does not imply a thermal spectrum. Then Tf is determined 
by the condition 

This shows that Tf ~ gV^ , at least with the approximation that we used for T pe ^ nl/ . A large energy density in relic 
gravitons gives a large contribution to the total density p and therefore to 17*. This results in a larger freeze-out 
temperature, more available neutrons and then in overproduction of 4 He. This is the idea behind the nucleosynthesis 
bound |73|. More precisely, since the density of 4 He increases also with the baryon to photon ratio 77, we could 
compensate an increase in with a decrease in 77, and therefore we also need a lower limit on rj, which is provided 
by the comparison with the abundance of deuterium and 3 He. 

Rather than g*, it is often used an 'effective number of neutrino species' N v defined as follows. In the standard 
model, at T ~ a few MeV, the active degrees of freedom are the photon, e ± , neutrinos and antineutrinos, and they 
have the same temperature, T, = T. Then, for N„ families of light neutrinos, g*(N u ) = 2 + (7/8) (4 + 27V„), where 
the factor of 2 comes from the two elicity states of the photon, 4 from e in the two elicity states, and 2N U counts 
the N v neutrinos and the N v antineutrinos, each with their single elicity state. For the Standard Model, N u — 3 and 
therefore g„ — 43/4. So we can define an 'effective number of neutrino species' N v from 

*wo = ? + E ^(jY+l E 4|) 4 - m 

i— extra bosons i— extra fcrmions 

One extra species of light neutrino, at the same temperature as the photons, would contribute one unit to N v , but 
all species, weighted with their energy density, contribute to N v , which of course in general is not an integer. For 
i = gravitons, we have gi = 2 and (Ti/T) 4 = p gw /Pj, where p 7 = 2(7r 2 /30)T 4 is the photon energy density. If 
gravitational waves give the only extra contribution to N v , compared to the standard model with N v — 3, using 
g*(N u ) = 2 + (7/8)(4 + 2N V ), the above equation gives immediately 
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, =l(N v -3), (47) 

H / NS 

where the subscript NS reminds that this equality holds at time of nucleosynthesis. If more extra species, not included 
in the standard model, contribute to g*(N v ), then the equal sign in the above equation is replaced by lower or equal 



The same happens if there is a contribution from any other form of energy present at the time of nucleosynthesis and 
not included in the energy density of radiation, like, e.g., primordial black holes. 

To obtain a bound on the energy density at the present time, we note that from the time of nucleosynthesis to the 
present time p gw scaled as 1/a 4 , while the photon temperature evolved following gs(T)T 3 a 3 = const, (i.e. constant 

entropy) and p 7 ~ T 4 ~ l/(a 4 <7g^ 3 ). Therefore 



Pgw \ / Pgw 



9s(T ) V /3 _/>gw\ /3.913\ 4/3 



Pi Jo V Pi /ns V3s(lMeV) ) \ p 7 J Ng V 10.75/ 
where the subscript zero denotes present time. Therefore we get the nucleosynthesis bound at the present time, 



(48) 



Pgw 

Pi 



< 0.227(AT, - 3) . (49) 



o 



Of course this bound holds only for GWs that were already produced at time of nucleosynthesis (T ~ MeV, t ~ 
sec). It does not apply to any stochastic background produced later, like backgrounds of astrophysical origin. Note 
that this is a bound on the total energy density in gravitational waves, integrated over all frequencies. Writing 
Pgw — / d(logf) dpgw/dlogf, multiplying both p gw and p 7 in eq. ( fl9| ) by h\j p c and inserting the numerical value 
h%p 7 /p c ~ 2.481 x 1CT 5 @, we get 

pf=oo 

/ d(bg/) h 2 n gw {f) <5.6x 10- 6 (iV,-3). (50) 
Jf=o 

The bound on N v from nucleosynthesis is subject to various systematic errors in the analysis, which have to do mainly 
with the issues of how much of the observed 4 He abundance is of primordial origin, and of the nuclear processing of 
3 He in stars, and as a consequence over the last five years have been quoted limits on N v ranging from 3.04 to around 
5. The situation has been recently reviewed in ref. [{741. The conclusions of ref. j74| is that, until current astrophysical 
uncertainties are clarified, N v < 4 is a conservative limit. Using extreme assumptions, a meaningful limit N v < 5 still 
exists, showing the robustness of the argument. Correspondingly, the right-hand side of eq. ( |50| ) is, conservatively, of 
order 5 x 10~ 6 and anyway cannot exceed 10~ 5 . 

If the integral cannot exceed these values, also its positive definite integrand /ig^gw(/) cannot exceed it over an 
appreciable interval of frequencies, A log / ~ 1. One might still have, in principle, a very narrow peak in /iQf2 gw (/) at 
some frequency /, with a peak value larger than say 10 -5 , while still its contribution to the integral could be small 
enough. But, apart from the fact that such a behaviour seems rather implausible, or at least is not suggested by any 
cosmological mechanism, it would also be probably of little help in the detection at broadband detectors like Virgo, 
because even if we gain in the height of the signal we lose because of the reduction of the useful frequency band A/, 
see eq. @). 

These numbers therefore give a first idea of what can be considered an interesting detection level for /igJl gw (/), 



which should be at least a few times 10 -6 , especially considering that the bound (50) refers not only to gravitational 
waves, but to all possible sources of energy which have not been included, like particles beyond the standard model, 
primordial black holes, etc. 



B. The astrophysical background 



GWs of astrophysical origin are not subject to the nucleosynthesis bound, and therefore our first concern is whether a 

gw( 



stochastic background due to a large number of unresolved astrophysical sources can give a contribution to /ip^gw(/) 



***To compare with eq. (56) of ref. |ll), note that while we use p 7 , in ref. the bound is written in terms of the total 
energy density in radiation at time of nucleosynthesis, which includes also the contribution of e , neutrinos and antineutrinos, 
(Prad)NS = [1 + (7/8) (2 + 3)](p 7 ) NS . 
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larger than the bound (p0|), or anyway larger than the expected relic signal, thereby masking the background of 
cosmological origin. 

A first observation is that there is a maximum frequency at which astrophysical sources can radiate. This comes 
from the fact that a source of mass M, even if very compact, will be at least as large as its gravitational radius 2GM, 
the bound being saturated by black holes. Even if its surface were rotating at the speed of light, its rotation period 
would be at least AttGM, and the source cannot emit waves with a period much shorter than that. Therefore we have 
a maximum frequency J75| , 

f<— 10 4 ^°Hz. (51) 

To emit near this maximum frequency an object must presumably have a mass of the order of the Chandrasekhar 
limit ~ I.2M0, which gives a maximum frequency of order 10 kHz |75f| , and this limit can be saturated only by very 
compact objects (see ref. (7(J] for a recent review of GWs emitted in the gravitational collapse to black holes, with 
typical frequencies /^S5 kHz). The same numbers, apart from factors of order one, can be obtained using the fact 
that for a self-gravitating Newtonian system with density p, radius R and mass M = p(4/3)7ri? 3 , there is a natural 
dynamical frequency [f77f 



/d^ = i(«Gp) 1/a =f™| (■">->• 



, 1/2 

With R > 2GM we recover the same order of magnitude estimate apart from a factor (3/8) 1//2 ~ 0.6. This is already 
an encouraging result, because it shows that the natural frequency domains of cosmological and astrophysical sources 
can be very different. We have seen that the natural frequency scale for Planckian physics is the GHz, while no 
astrophysical objects can emit above, say, 6-10 kHz. A stochastic background detected above these frequencies would 
be unambiguosly of cosmological origin. 

However, ground based interferometers have their maximum sensitivity around 100 Hz, where astrophysical sources 
hopefully produce interesting radiation (since these sources were the original motivation for the construction of 
interferometers) . The radiation from a single source is not a problem, since it is easily distinguished from a stochastic 
background. The problem arises if there are many unresolved sources. (Of course this is a problem from the point of 
view of the cosmological background, but the observation of the astrophysical background would be very interesting 
in itself; techniques for the detection of this background with a single interferometer using the fact that it is not 
isotropic and exploiting the sideral modulation of the signal have been discussed in ref. ]7q|). 



The stochastic background from rotating neutron stars has been discussed in 1 79 and references therein. The main 
uncertainty comes from the estimate of the typical ellipticy e of the neutron star, which measures its deviation from 
sphericity. An upper bound on e can be obtained assuming that the observed slowing down of the period of known 
pulsars is interely due to the emission of gravitational radiation. This is almost certainly a gross overestimate, since 
most of the spin down is probably due to electromagnetic losses, at least for Crab-like pulsars. With realistic estimates 
for e, ref. ]79| gives, at / = 100 Hz, a value of h c (f) ~5x 10~ 28 , that, using eq. (^9|), corresponds to /igft gw (100Hz) ~ 
10 -15 . This is very far from the sensitivity of even the advanced experiments. An absolute upper bound can be 
obtained assuming that the spin down is due only to gravitational losses, and this gives /igO gw (100Hz) ~ 10~ 7 , but 
again this value is probably a gross overestimate. 

The stochastic background from supernovae is studied in ref. [§0j . The expected frequencies in this case are of the 
order of the kHz or lower (down to 500-600 Hz, depending on the redshift when these objects are produced). Using the 
observational data on the star formation rate, it turns out that the duty cycle, i.e. the ratio between the duration of a 
typical burst and the typical time interval between successive bursts is low, of order 0.01. Therefore this background 
is not stochastic, but rather like a 'pop noise', and can be distinguished from a really stochastic background. The 
value of /ipfigw for the background from supernovae have been computed in [{sof assuming axially symmetric collapse, 
and assuming that all sources have the same value of a — J/(GM 2 ), where J is the angular momentum. The results 
depend on the value of a, and on ho. For typical choices, one gets values of order /igf2 gw ~ 10~ 14 — 10~ n at say 
/ = 65 Hz, rising up to fogfigw ~ 10~ 9 — 10~ 8 around 3kHz, where the spectrum has the maximum [ |S0"| |. 

These results suggest that astrophysical backgrounds might not be a problem for the detection of a relic background 
at LIGO/ Virgo frequencies. The situation is different in the LISA frequency band ]7^ , [77| , ^l| , [79|Jl9| ] . LISA can reach a 
sensitivity of order ft.gf2g W ~ a few xlO -13 at / ~ 10 -3 Hz (see fig. 1.3 of ref. |lj|). However, for frequencies below a 
few mHz, one expects a stochastic background due to a large number of galactic white dwarf binaries. The estimate 
of this background depends on the rate of white dwarf mergers, which is uncertain. With rates of order 4 x 10 -3 per 
year (which should be a secure upper limit [p9[), the background can be as large as /ig£l gw ~ 10~ 8 at / = 10~ 3 Hz. 
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This number is actually quite uncertain, and in ref. |^9|, fig. 1.3, it is used another plausible rate, which gives for 
instance /ij^gw ~ 10~ n — 10 -10 at / = 10~ 3 Hz. Above a frequency of order a few times 10~ 2 Hz, most signals from 
galactic binaries can be resolved individually and no continuous background of galactic origin is presently known at 
the level of sensitivity of LISA. 

It should be observed that, even if an astrophysical background is present, and masks a relic background, not all 
hopes are lost. If we understand well enough the astrophysical background, we can subtract it, and the relic background 
would still be observable if it is much larger than the uncertainty that we have on the astrophysical background. In 
fact, LISA should be able to subtract the background due to white dwarf binaries, since there is a large number of 
binaries close enough to be individually resolvable jlj|. This should allow to predict with some accuracy the space 
density of white dwarf binaries in other parts of the Galaxy, and therefore to compute the stochastic background that 
they produce. Furthermore, any background of galactic origin is likely to be concentrated near the galactic plane, and 
this is another handle for its identification and subtraction. The situation is more uncertain for the contribution of 
extragalactic binaries, wheih again can be relevant at LISA frequencies. The uncertainty in the merging rate is such 
that it cannot be predicted reliably, but it is believed to be lower than the galactic background Q . In this case the 
only handle for the subtraction would be the form of the spectrum. In fact, even if the strength is quite uncertain, 
the form of the spectrum may be known quite well [fl9| . 



C. Cosmological predictions for the intensity 

The next question is whether it is reasonable to expect that some cosmological production mechanism saturates the 
nucleosynthesis bound. This of course depends on the production mechanism, but some general considerations are 
possible. First of all, it is clear that, if GWs are produced at the Planck scale by collisions and decays, togheter with 
the photons that we observe today in the CMBR, and there is not an inflationary phase at later time (the scenario of 
sect. 6B), we expect roughly p gw ~ p 7 , so that the bound (E^) is approximately saturated. 

If the mechanism that produces GWs is different from the mechanism that produces the photons in the CMBR, 
often it is still possible to relate the respective energy densities, simply because the scales of the two processes are 
related. For instance, ref. |36[ discusses the spectrum of relic GWs produced in string cosmology, using only general 
arguments. The peak of the spectrum, fx, is fixed with the criterium iJ» ~ Ms, which as discussed in sect. 5, 
corresponds to an effective 'temperature' at time of production T* ~ (M^Mpi) 1 ' 2 = (Ms/Mpi) 1 / 2 Mpi. Redshifting 
T* we get a characteristic frequency today fx ~ {Ms/My>\) 1 / 2 Tq. Fixing the peak energy from the 'one graviton 
level' Q, iif ~ 1, eq. ( |2l| ) gives a peak energy density dp/d\ogf ~ ff, so that at the peak frequency we have 

dp g w ( M s \ dp y 

dlogf ~ \M P J dlog/' 1 ) 

so that p gw is related to p 1 , although with a suppression factor (M$ /Mp\) 2 which is expected to range between 10~ 4 
and 10~ 2 , and numerical factors. 

Similarly, in the production of GWs through bubble collisions when inflation terminates with a first order phase 
transition, there is only one energy scale, the vacuum energy density M during inflation (so that M 4 is the false- 
vacuum energy density). The value of M fixes the reheating temperature, and therefore p 7 , from (ir 2 <u/30)T^ h ~ M , 
and fixes also the energy liberated in gravitational waves. The latter can be estimated as follows |6lj. The typical 
wavelength A* produced in the collision of two bubbles will be of the order of the radius of the bubbles when they 
collide, which in turn is a fraction of the horizon scale H' 1 . As in sect. 4, we write A* = eif" 1 , and ~ M 2 /Mp\. 
The energy liberated in GWs in the collision of two bubbles is of order Eqw ~ GM|/A» where Mb ~ M 4 A 3 is the 
energy of a typical bubble. The fraction of the false vacuum energy that goes into GWs instead of going into p 7 is 
therefore Eqw/Mb ~ e 2 - This gives, at time of production 

p g w ~ <?p 1 ■ (54) 

Again p gw is related to p 1 , and the basic reason is that there is essentially only one dimensionful parameter that 
enters in the estimate. 

From these explicit examples, we see that independently of the production mechanism, eq. ( |54| ) is quite general. 
The scale for the intensity of a relic GW background is indeed fixed in many cases by p 7 , which therefore gives 
a first order of magnitude estimate of the effect. However, there are also suppression factors, like the factor e 2 in 
eq. ( ]54| ) or (Mg/Mpi) 2 in eq. (|53|), that, togheter with the exact numerical coefficients, are crucial for the detection at 
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present experiments. A value e ~ 0.1 would allow detection at the level /ig^gw ~ 10~ 7 while e ~ 10~ 2 would require 
/ipfigw ~ 10~ 9 , which is beyond the possibilities of first generation experiments. 

For GWs produced by collisions at the Planck era we can even expect e ~ 1. However, another type of suppression 
is present if we have a spectrum of relic GWs with a total energy density p gw ~ p 7 with a cutoff in the GHz, and 
we want to observe it in the kHz region. In this case, as we discussed, we must hope that the spectrum is practically 
flat between the kHz and the GHz. While we have seen in sect. 5 that there are good reasons for considering such 
a behaviour, the fact that the energy density is spread over a wide interval of frequencies diminuishcs of course the 
value of h^fl gw (f) at a given frequency. The nucleosynthesis bound ( |50| ) then gives a maximum value at 1kHz 

frfo gw (lkHz) < 5 " 6 * ^i^r 3) - 4 x l(T 7 (iV„ - 3) . (55) 

If the spectrum extends further toward lower frequencies, the maximum value of h^Clgv, decreases accordingly, with a 
factor log(lGHz// m in) instead of log(lGHz/lkHz). 

In production mechanisms where there is no obvious relation between p gw and p 7 , there is no general reason 
suggesting a value of p gw close to the nucleosynthesis bound, but still we can have rather large relic backgrounds. An 
important example is provided by cosmic strings (see |^| and references therein). Cosmic strings might form during 
phase transitions, but whether they actually form or not depends on the precise nature of the transition. Therefore, 
as remarked in ref. [ [TT|| , the computation of the GW spectra from cosmic strings must be considered as illustrative 
rather than realistic (as on the other hand happens in most of the examples that we have discussed). Anyway, at 
LIGO/Virgo frequencies, the typical estimate of the intensity is of order ti$l gw ~ 10~ 8 — 10~ 7 . (Similar results are 
obtained from hybrid topological defects |^| ). The scale for these numbers is given first of all by the combination 
(Gp,) 2 , where p is the mass per unit length of the string. It comes from the fact that a loop radiates with a power 
P ~ "fGp 2 , (7 is a dimensionless constant) while, evaluating ft,§£l gw , another factor of G comes from l/p c - For strings 
created at the GUT scale Gp ~ 10 -6 . Then, there are also large numbers related to the number of strings per horizon 
volume, which is thought to be of order 50, to the fact that also 7 ~ 50, and to the size of the loop at formation 
time pH , that finally give a value /igf2 gw ~ 10~ 8 — 10~ 7 . It is clear that in this case the right order of magnitude can 
only be obtained with a detailed understanding of the dynamics of the cosmic string network. 



VIII. CONCLUSIONS 



Present GW experiments have not been designed especially for the detection of GW backgrounds of cosmological 
origin. Nevertheless, there are chances that in their frequency window there might be a cosmological signal. The most 
naive estimate of the frequency range for signals from the very early Universe singles out the GHz region, very far 
from the region accessible to ground based interferometers, / < a few kHz, or to resonant masses. To have a signal in 
the accessible region, one of these two conditons should be met: either we find a spectrum with a long low-frequency 
tail, that extends from the GHz down to the kHz region, or we have some explosive production mechanism much 
below the Planck scale. As we have discussed, both situations seem to be not at all unusual, at least in the examples 
that have been worked out to date. The crucial point is the value of the intensity of the background. 

With a very optimistic attitude, one could hope for a signal, present just in the 10Hz-l kHz band, with the maximum 
intensity compatible with the nucleosynthesis bound, /iQf2 gw ~ a few xlO -6 (or even 10~ 5 , stretching all parameters 
to the maximum limit). Such an option is not excluded, and the fact that such a background is not predicted by 
the mechanisms that have been investigated to date is probably not a very strong objection, given our theoretical 
ignorance of physics at the Planck scale and the rate at which new production mechanisms have been proposed in 
recent years, see the reference list. However, with more realistic estimates, on general grounds it appears difficult to 
predict a background that in the kHz region exceeds the level h^fl gw ~ a few xl0~ 7 , independently of the production 
mechanism. This should be considered the minimum detection level at which a significant search can start. Such a 
level is beyond the sensitivity of first generation experiments, unless a ground based interferometer is correlated with 
a second interferometer, located at a distance small enough so that a significant correlation is possible, and large 
enough to decorrelate local noises. A few tens of kilometers would probably be the right order of magnitude. In this 
case we could detect a signal at the level 2 x 10 -7 in one year of integration time, with SNR=1.65, and the level of a 
few x 10 -8 could be reached with longer integration time (and possibly allowing for a slightly worse confidence level; 
this could make some sense in a stochastic search because the SNR increases with time, and the hint of a signal at 
low SNR would provide a strong motivation for pursuing the search). The difficulties of such a detection are clear, 
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but it should be stressed that the payoff of a positive result would be enormous, opening up a window in the Universe 
and in fundamental high-energy physics that will never be reached with particle physics experiments. 
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APPENDIX A: 

In this Appendix we discuss the estimate of the quantity e, that enters eq. d36|), in two specific examples. The 
discussion is meant to illustrate the uncertainty inherent in the determination of e and the range of values that it can 
take. 

i) Inflation- RD transition. Quantum particle creation due to vacuum fluctuations is associated with a sudden 
change of the scale factor, e.g. from a DeSitter inflationary phase to the RD phase Ji^-^(|. Let us consider a typical 
inflationary scenario [H: at some time A* the potential energy V(<j>) of some scalar field (p (the 'inflaton' field) starts to 
trigger inflation, and the Universe expands exponentially. This inflationary stage ends at time tf, when the so-called 
slow- roll conditions on V{<f) are not anymore satisfied. To solve the problems of standard cosmology, tf — U should 
be at least of order 60i/pg, where Ads is the Hubble constant during the DeSitter inflationary phase. At t = tf the 
temperature of the Universe is exponentially small, and a reheating period must take place. In the simplest scenarios, 
reheating lasts for a time interval At ^ T^ 1 , where is the decay width of the inflaton, and the decay of the inflaton 
reheats the Universe, which then enters the standard RD phase. 

The cutoff of the spectrum at large frequency is determined by how fast the transition between the two phases takes 
place |4[],|^|. Therefore the maximum value of the frequency produced is /* ~ 1/At ~ T^. The gravitons have been 
produced at a time i*, with tf < i* < tf + At. If the transition proceeds slowly, At ^> tf, i.e. if is sufficiently 
small compared to ~ H^s, then tf + At ~ At and we can take i* ~ At ~ r^ 1 . Therefore in this case the statement 
/* ~ H(t*) = l/(2t») is indeed correct. However, if is sufficiently large, so that the transition takes place on an 
interval At <g; tf, the value of time explored detecting these gravitons is t* ~ tf, while the characteristic frequency 
produced is still /* ~ T^. In this case /* ^> 1/f*, or A* <C H^ 1 . 

However, in most models <C -Hds- F° r instance, in models derived from supergravity, a large would produce 
a large reheating temperature T r h, and a value T r h > 10 9 GeV is forbidden because of the gravitino problem (see e.g. 
ref. ||, pg. 298). And indeed small values of come out quite naturally in supersymmetric models. For instance, 
in the model of ref. Q, ~ A 6 /M 5 where M — A/pi/\/87r and A~3x 10~ 5 M is a mass scale presumably related 
to supersymmetry breaking. Therefore in this case ~ 10~ 27 M is extremely small. 

Not all models of inflation satisfy <C i?DS automatically. For instance, in the old Coleman- Weinberg model 
> -HdSj but the model is not viable. So, in the case of quantum production of particles after an inflationary phase, 
the estimate /* ~ 1/t*, and therefore e ~ 1, is correct, at least in typical viable models. 

it) Phase transitions. Another important example is given by phase transitions in the early Universe, as the 
QCD or the electroweak phase transition. These are believed to proceed through bubbles nucleation of the low- 
temperature phase, when the Universe cools below the critical temperature for the phase transition, T c . If the 
transition occurs explosively, the collisions of different bubbles produce gravitational waves @,|5) ■ The characteristic 
wavelength of these GWs is of the order of the typical radius R of the bubbles when they collide. We write again 
Apeak it*) ~ R = tH^ 1 . Plausible values of e have been discussed in detail by Hogan and by Witten (ref. [jS5| , 
app. B). In the case of the QCD transition e ~ 1 is excluded because otherwise in the bubble collisions there would 
be overproduction of primordial black holes. Production of primordial black holes is severely constrained because 
the energy density of primordial black holes scales like 1/a 3 in the RD phase, and therefore they would come to 
dominate the energy density of the Universe. In particular, the Universe would not be radiation dominated during 
nucleosynthesis. Actually, primordial black holes can evaporate through emission of Hawking radiation. However, in 
order to evaporate before nucleosynthesis, the scale at which the phase transition takes place must be higher than 
10 11 GeV Therefore e could be of order one in phase transitions which take place at these temperatures, but not 
in the QCD or electroweak phase transition. 

To estimate e it is necessary to make assumptions about how the phase transition is nucleated. If it is nucleated 
by thermal fluctuations, ref. |84j suggests an upper bound e<10 -2 . For nucleation of bubbles via quantum tunneling 
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a detailed analysis has been done in ref. [p6| . The relevant parameter is the nucleation rate, which again is a model- 
dependent quantity, but for a wide class of models one typically finds e ~ 10 -3 — 10~ 2 ]8rj , |87| , |3l| ] . If the transition is 
nucleated by impurities (which is most often the case, except in very pure and homogeneous samples) the issue is much 
more complicated. For instance, the impurities could be given by turbolent motion of the cosmic fluid generated prior 
to the QCD epoch; they would depend on the detailed spectrum of the fluid motion, or in general on the characteristic 
distance between impurities, and in this case e is basically impossible to estimate, even as an order of magnitude. 
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